Chromoelectric response functions for quark-gluon plasma 
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We determine the chromoelectric response of quark-gluon plasma (QGP) systematically within 
the framework of classical transport equations. The transport equations are set up in the phase 
space which includes the SU (3) group space corresponding to color (which is a dynamical degree of 
freedom), in addition to the position - momentum variables. The distribution functions are defined 
by projecting the density operators for the quarks and the gluons to their respective coherent 
states (defined over the extended phase space). The full import of the Yang-Mills(YM) dynamics 
is shown to manifest through the emergence of an intrinsic nonlinear, nonlocal response, whose 
behavior we determine in the long wavelength limit. It also manifests as a tensor response which is 
a characteristic of gluons. The response functions are shown to have a natural interpretation in terms 
of the renormalizations of the Abelian and the non-Abelian coupling constants. A detailed analysis 
of the screening of heavy quark potential and of the exact role played by the Debye mass screening 
in the case of the Cornell potential, is performed. We also discuss the non-Abelian contribution to 
Landau damping in QGP. 
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I. INTRODUCTION 



The purpose of this paper is to study systematically the 
color response functions of quark gluon plasma, with a 
proper incorporation of the non-Abelian dynamics. Con- 
sidering the color electric case, we show the emergence of 
additional response functions which do not have counter- 
parts in the well studied electrodynamic plasmas. More 
precisely, we show the emergence of a non-Abelian com- 
ponent of the chromoelectric permittivity in the matter 
and gluonic sectors. The gluonic sector is shown to ex- 
hibit yet another response, corresponding to its color- 
tensor exciatations. In this work, we merely illustrate 
these responses in the simple case of ideal plasma. We 
pay particular attention to the role of screening length in 
modifying the dynamics of heavy quarkonium systems. 
Applications to more realistic cases will be taken up in 
a separate work. The analysis is performed within the 
framework of a classical transport equation - the anlog 
of Vlasov equation - as is appropriate to non-Abelian 
dynamics. We start with the motivation below. 

Searches for QGP in Ultra Relativistic Heavy Ion Col- 
lisions (URHIC) have been on experimentally [l|, 0, 0, Q 
for more than a decade. We have today a wealth 
of information coming from observations involving (i) 
particle multiplicities, (ii) minijcts, (iii) jet quenching, 
(iv) strangeness enhancement, (v) collective flows, (vi) 
heavy quark dynamics (J/^* suppression) Q , (vii) Han- 
bury Brown- Twiss interferometric measurements |6| etc. 
It is strongly suspected that QGP is already produced 
in experiments at RHIC. Further, it has been inferred 
from flow measurements that the deconfined phase is, in 
all likelihood, a liquid state with a very low viscosity - 
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smaller than even that of liquid helium 0, 0] . The up- 
coming experiments at Large Hadron Collider(LHC) may 
be expected to reveal many more (unexpected) results. 

A parallel development that has taken place theoret- 
ically is the realization that QGP exhibits a strong col- 
lective behavior which cannot be understood perturba- 
tively [1, Q. First of all, improved lattice computations 
[9l.[l0| predict that the deconfined phase close to the tran- 
sition temperature Tc is highly non-perturbative, with 
the equation of state being far away from that of the 
ideal gas of quarks and gluons. The ideal behavior is not 
expected to be realized even at T = 2Tc. Analytic stud- 
ies [iJI based on improved perturbative approaches (with 
computations of 0(.g^log(g)) also arrive at the same con- 
clusion and indicate that the ideal behavior will be seen 
only at much higher temperatures. Experimentally, the 
flow measurements at RHIC do indicate a low viscosity 
liquid like behavior of the QGP [T3| • Finally, it is known 
by now [l^, [3, [iBl that a classical behavior emerges nat- 
urally when one considers hard thermal loop(HTL) con- 
tributions. A local formulation of HTL effective action 
has been obtained by Blaizot and lancu who have suc- 
ceeded in rewriting the HTL effective theory as a kinetic 
theory with a Vlasov term [H, [l^, [Hi • 

It may be, therefore, a fruitful endeavor to take the 
classical frame work seriously and explore the extent to 
which it can capture the expected features of QGP. It 
is well recognized that its success would depend cru- 
cially on the availability of reasonable mechanisms for 
(i) the production of soft partons (rate term in the phase 
space) and (ii) the expression for the equation of state 
and its equilibrium distribution function. Indeed, sig- 
nificant progress has been made in both the directions 
by employing a variety of techniques [H, [13, [Ml, [ll] , al- 
though they are all not always mutually consistent. Be 
it as may, we wish to point out in this paper that the 
evaluation of the response functions of the QGP, so cru- 
cial for the signals, is further dependent on the form of 
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the distribution functions in the color space, which is in- 
herited from the parent density operators; it may not be 
assigned at will. 

The central theme of the paper is to couple the coher- 
ent state representation in the color space with the YM 
dynamics in the transport equations. In doing so, we de- 
rive, in a consistent manner, the response functions which 
show manifestly the Abelian and the non-Abelian contri- 
butions. The latter is nonlinear and nonlocal in space 
time. Hence it is also non-Markovian. We also show the 
intrinsic distinction between a pure QCD plasma which is 
purely gluonic and QGP which has both matter and glu- 
onic components. Since we do not wish to make any real- 
istic calculation here, many simplifications will be made 
in the transport equations: (i) the equilibrium distribu- 
tion functions will be taken to be ideal, (ii) we study the 
response when the system is close to equilibrium, so that 
the source term may also be dropped and (iii) finally, we 
take the plasma to be spatially isotropic. The plasma 
is taken to be collisionless. These simplifications will be 
improved upon in a separate paper, where more realistic 
collision terms such as the one derived by Bodeker [2^ 
will be used. As remarked above, the primary purpose 
of the paper is to elucidate the structure of the response 
functions and to illustrate it in simple situations. 

The chromoclectric response functions are particularly 
important, ever since Matsui and Satz predicted J/^ 
suppression as a signature for QGP in heavy ion colli- 
sions. Lattice calculations @ do support the prediction; 
early estimates of the suppression in heavy ion collisions 
have involved the assumption of a hydrodynamic expan- 
sion of an ideal plasma. The suppression is of course 
determined by the Debye mass which evolves with the 
temperature. We study the precise role of the Debye 
mass in the specific case of the Cornell potential. In any 
case, it appears that the Debye mass is an incomplete 
manifestation of the non-Abelian dynamics, in as much 
as that it does not require the non-Abelian interactions 
for its emergence. It is not unreasonable to look for exclu- 
sively non-Abelian responses which would be nonlinear, 
and reflect additional properties of chromodynamics in a 
medium. We do show that such a non-Abelian permittiv- 
ity exists. The new permittivity may also lead to other 
signatures which can perhaps be tested experimentally. 



A. A brief review 

We briefly review the studies on the QGP response 
functions so far to the extent that they are relevant to 
our work. There have been many attempts to calcu- 
late the response function for QGP. An approach based 
on perturbative QCD (pQCD) at finite temperature was 
adopted by Weldon [2^ who studied plasma screening 
and plasma oscillations. Based on an Abelian analysis, 
Mustafa, Thoma and Chakraborty [2^ have argued for 
bound quark structures, which have also been predicted 
by Liao and Shuryak [23] ■ The conclusion is essentially 



based on their analysis for moving partons. Petreczky 
[2^ has, on the other hand, considered a SU{2) plasma, 
and extracted the chromoclectric screening by analyzing 
the long distance behavior of the static part of the lon- 
gitudinal propagator in the lowest order. His results are 
off by about 25% from lattice simulations, and are not 
easily adapted to the scenario prevalent in URHIC. Fi- 
nally, in a work somewhat close to ours, Chen et al poj 
have employed the Vlasov equation to study the color re- 
sponse. However, the gauge invariance of their results is 
not apparent, as also its applicability to URHIC. None of 
the above approaches employs the coherent state projec- 
tions which establish a natural connection between the 
underlying quantum dynamics with the classical approxi- 
mations. They do not also study exhaustively all possible 
response functions involving the color degree of freedom. 

II. DETERMINATION OF THE RESPONSE 
FUNCTIONS 

As mentioned above, the approach will be based on 
the transport equation for the phase space distribution 
functions for the quarks and the gluons. The distribution 
functions for an N particle system will have the generic 
form /(qi, • • • qXr; pi,---pN; Q?, •• -Q?/) in terms of the 
coordinates q, the momenta p, and the color charges 
which are defined in the space corresponding to the gauge 
group. Note that the distribution function is invariant 
under gauge transformations (for a formal demonstra- 
tion, see e.g., (30l|). Let us first consider the generic form 
of the transport equation 

W + + + ^ 90^ = ^ + ^' 

where the source term is denoted by E and C is the colli- 
sion term. The source term is indeed important in study- 
ing the production and evolution of QGP [3l|, [H, l3a |. 
But as mentioned above, we study the response when 
system is close to its equilibrium configuration, whence 
the source term may be dropped. We further work in 
the Vlasov limit so that the collision term will also be 
dropped. The responses are, of course eventually evalu- 
ated in the limit that the collision time Tc 00, via the 
standard Landau prescription. The reason for working 
in this limit is that the inclusion of any realistic colli- 
sion term, as, for instance, the one derived by Bodeker 
[2^ will complicate the discussions. On the other hand, 
a simple minded relaxation term would violate current 
conservation [3^ . This draw back will be remedied in a 
separate paper. 

The third term in EqlT] represents the Vlasov term 
corresponding to the action of the mean field on the 
plasma. The next term, involving the derivative of the 
color charge, is unique to the non-Abelian plasma, and 
gives the dynamical variation of the color charge (in 
the compact part of the phase space). The dynamics 
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modulates and modifies the contribution of tlie standard 
Vlasov term; furtlier, it also gives rise to the non-Abelian 
component of the response function. Our analysis makes 
crucial use of the most general form of the quark and the 
gluonic distribution functions in the color variables, and 
we take up the determination of their form in the next 
section. 



A. The distribution functions 

1. The coherent bases 

We elucidate the procedure for obtaining the generic 
forms of phase space distribution functions from their 
underlying quantum states. The best way of extracting 
them is to project the parent quantum state to a coher- 
ent basis. Indeed, with the usual position-momentum 
variables, the coherent states are the closest to the clas- 
sical states since they possess minimum uncertainty in 
position and momentum. The projection defines, in a 
natural manner, the area over which the phase space is 
coarse grained: it is simply given by 27r?i for each de- 
gree of freedom. Further, the coherent state projection is 
equivalent to the Wigner distribution function used e.g. 
by Elze and Heinz |35l. [36t . However, it has the advantage 
that unlike the Wigner distribution function which can 
take negative values, the phase space density obtained 
from the coherent state is always non-negative. 

The above method needs a refinement when it comes 
to the color degree of freedom. Recall that the classical 
phase spaces associated with finite dimensional Hilbert 
spaces are compact (as for example, in the case of spin) . 
There is no neat separation of canonically conjugate vari- 
ables, and there is no associated commutator (of the kind 
[x,p\ = i) that leads to minimum uncertainty states. In- 
stead, the commutators are given by the Lie Algebra of 
the associated compact group. Care needs to be taken 
in defining the phase space distribution in this case. The 
most convenient method is to follow Perelomov [33] : take 
a representative state, and act the group exp{iTa9a} on 
a reference state which is typically taken to be the state 
with the highest weight. The resulting orbit, 0(|{6'a} > ) 
forms a faithful copy of the associated phase space. The 
quantity {9i, ■ ■ ■ 9n\p\0i, ■ ■ ■ On) is interpreted consistently 
as the classical probability density in the color part od 
the phase space. These distributions so obtained are al- 
ways smeared, for any finite dimensional representation 
of the compact group. 

2. Distribution functions for quarks and gluons 

We employ the method outlined above for the specific 
system of interest. Let p be the density operator for a 
parton. Denote by j^* >c= \r, p> (8)|(5° > the coherent 
basis for the parton, where the first term is the usual 
minimum uncertainty state. As mentioned, the latter 



term is obtained by the action of the gauge group on 
a standard state, \ip >s, taken to be the one with the 
highest weight [33| . Thus the coherent states in the color 
sector have the form jQ" >= exp{iQ'^T°'}\ips >, where 
the variables provide a coordinate description of the 
group space. For example, the parameters for SU{2) may 
be chosen to be the Euler angles. A generalized Euler 
angle description for SU{3), which is of relevance to us 
here has been provided by Byrd [111 . We will be using it 
subsequently. 

The form of the generators are to be chosen de- 
pending on the parton, i.e., the representation to which 
it belongs. The quarks transform according to the funda- 
mental representation of the gauge group and the genera- 
tors will also be chosen accordingly. Gluons, on the other 
hand, belong to the adjoint representation. As a direct 
consequence, the dimension of the phase space will be 
six for the quarks, and eight for the gluons if we consider 
SU{3) gauge group. Incidentally, note that the phase 
space for the quark is isomorphic to its Hilbert space. 
The gluonic phase space constitutes, in contrast, a coset 
space, being an orbit of the gauge group in the Hilbert 
space which is of dimension fifteen. The property of the 
coherent state depends naturally on the orbit to which 
it belongs. Suffice to say, the color coherent states are 
natural generalizations of the more well known s pin co- 
herent states. We refer the reader to literature [33| for 
more details. 

The structure of the copy of the phase space, obtained 
thus as a coset space has strong implications on the na- 
ture of the color distributions, bearing a deep connection 
with the classical-quantum correspondence. For any fi- 
nite dimensional representation of the group, only a finite 
number of the moments of the variables can contribute. 
To illustrate with a simple example, the state of a spin 
half particle is characterized entirely by the vector polar- 
ization Tt{pS}. a spin one state requires, in addition, 
a specification of the second order tensor polarization 
Tt{pS,,}, where 5„ = ^i^i+Mi _ ^s^S.j. 

We perform a similar analysis for quarks and gluons. 
Since the quarks belong to the fundamental represen- 
tation of the gauge group, the quark distributions can 
at the most be linear in the color charge Q. All the 
higher order terms vanish identically. Thus, the single 
quark distribution function has the most general form 
/ = /o + /"O" where /o, are functions of the usual 
phase space variables {xi, pi}. In contrast, the gluons 
belong to the adjoint representation; they admit a more 
general expansion, with an additional bilinear term in the 
color charge. Note that the Abelian limit is obtained by 
dropping all the multipole terms except the scalar com- 
ponent. The classical limit is obtained in the other limit, 
the so called large Nc limit, by admitting multipole con- 
tributions of arbitrary high orders. Finally, it may be 
noted that the two Casimirs Q^Q" and d"^" Q" Q" re- 
strict the dynamics under the gauge independent inter- 
actions to a smaller subspace. As we shall see latter, the 
above resolution of the distribution functions brings out 
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the richness m the response functions. 



III. THE RESPONSE FUNCTIONS 

A. The quark-ant iquark sector 

Strictly speaking, one should write a two particle ditri- 
bution function for the qq system to study the color 
charge excitations. This procedure is laborious (39j and 
we present a simplified version below. Consider color 
deviation of the qq system from its equilibrium config- 
uration /o, which is taken to be a uniform distribution 
in color space. Accordingly, we write the distribution 
function in the form 



(2) 



which emphasizes that /o is independent of position and 
color variables. The color fiuctuations which are obtained 
by perturbing around /o are denoted by Q°'f°'- The color 
density functions are functions of both position and 
momentum. Finally, we also note that at equilibrium, 
the chromoelectric field vanishes identically. 
The transport equation reads 



dt ^ dxi ^ dpi dQ"- 



0, 



which on the application of the Wong equation [40 



dT 



(3) 



(4) 



acquires the form 



-/'-"{A[,(r)-z;.A^(f)}Q'" 



Of 



0, 



(5) 



in obtaining which we have assumed that the plasma is 
spatially isotropic and that the chromomagnetic field is 
absent. Unlike the electrodynamic plasma, the above 
equaton needs a specification of the gauge 41 1 [i^ . We 
study EqlHin the temporal gauge, Aq = 0. In this gauge 
the chromoelectric field has a simple dependence on the 
gauge potential: 



dt 



The transport equation now gets a simpler form 
9f . 



5/ 
dt 



' dxi 



The evaluation of the chromoelectric response requires 
the determination of the equations for the moments of 
the distribution function, involving integrations in the 



group space by employing the appropriate Haar measure 
in the group space. The details of the evaluation of the 
color moments are provided in the Appendix. 

At this juncture, it is convenient to switch over to the 
vector notation for color variables: 

/°Q° = / • Q; = (i* X By\ 

The latter notation collapses to the usual notation when 
the gauge group is SU{2). We continue to denote the 
position and momenta by their components. Thus, we 
write 



,A, X / = 0. 



(7) 



Note that in obtaining the above equation we have fur- 
ther assumed that / is a small perturbation on /q. In 
particular, dominates over the gradient of the color 

vector fluctuation |^| which we drop. 

The color response function is to be derived from Eql7| 
which, as we see, simply replaces the ordinary deriva- 
tive by the covariant derivative, in the temporal gauge. 
It turns out that the extraction of the response func- 
tion is not as straight forward as in the case of ordinary 
plasma. We follow the standard treatment and write 
Eql7]in the Fourier space. We denote by 4'{u)^ ki,Pi), the 
Fourier transform of f{t, Xi,pi). Similarly, £i and At rep- 
resent the Fourier transforms of the gauge field and the 
gauge potential respectively. We obtain 



{Lo,ki,pi) = ~Q£i-^^i^ + {Ai-k $){uj,ki,pi)^Vi. (8) 



V 



In the above expression, T) = iuj — ikiVi. We have em- 
ployed the notation ★ to denote the convolution of the 
Fourier transforms, combined with the cross product in 
the color variables: 

{Ai^Bf = r''' J duj'd?k'A\Lo - Lo', k, - K)B\uj', k[). 

Q denotes the magnitude of the color charge. 

The expression for the charge density p is easily ob- 
tained to be 

p{LO,k,) = -Q''S^ J ^'P^§^ + Q{A*^^){uJ,h), (9) 

where we have defined 

$i{u;,ki) = J d^pvi^${uj,k,,p^). (10) 

It is clear from the above equation that the expression 
for p is not closed, even after we employ the Gauss Law 
expression, all beacuse of the second term which cannot 
be written in terms of the charge density or the field. A 
rigorous evaluation would lead to an infinite hierarchy 
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of equations with moments of the distribution functions 
with respect to velocities. Since we are interested in the 
bulk properties of the medium, we evaluate Eq|H] in the 
long wavelength limit, by assuming that ki/ui <C 1. In 
this approximation, Eq. 1101 simplv becomes 

ijj 

in terms of the color current density J^. It follows from 
Eq.8 that the equation for Ji is given by 



where 



$,j{uj,ki) = j 



= I d pViVj—(l){uj,ki,p,), 



(11) 



which again involves an expression which is one order 
higher in the moment with respect to velocity. We trun- 
cate the hierarchy, by dropping the contribution from the 
second term, and write 



Q tj / d^pvt — ——. 



We may thus rewrite Eq|5] in the form 

p{u!, h,pi) = -iQ^EJi - Q^—Ai ★ 



(12) 



(13) 



The integrals are entirely functions of the equilib- 

rium distribution functions, and given by 



,3 15/o 

,3 1 dfo 
d PVi^-^ 



(14) 



The determination of permittivity is almost complete. 
Since we are considering an isotropic system, we observe 
that the integrals above have the resolutions 



h ^ kiIo{Lo,k) 

It] = Stjli{uj, k) + kikjhiuj, k). 



(15) 



Clearly, I2 is analytic in k and, thus does not contribute 
at low wavelengths. The final expression is thus given by 

p{uj, h) + iQ'^ki£t{u;, ki)k){uj, k) Ai * (Eih) = C(.16) 

UJ 

The response functions can be inferred by comparing 
Eql 131 with the Gauss' law in vacuum in the Fourier space: 



p{uj, ki) — iki£i{uj, ki) + Ai * £i — 0. 



(17) 



Combining EasfT51 and [T71 we are led to define the fol- 
lowing permittivities 

et'iiu, fc) = {1 + Q^IoiuJ, fc)}%<5''^ = eAStjSab. (18) 



e^^^iuj,uj') = {l + 



Q'hito\k')\,, 



(19) 



Note the emergence of the new permittivity ejv which 
has no Abelian counterpart. It is an invariant tensor of 
rank three in the color space (being proportional to f^^^). 
The new permittivity emerges, we stress, in addition to 
the Abelian permittivity e^, which leads to the standard 
screening. As much as the Abelian permittivity signifies 
a renormalization of the charge, the non-Abelian per- 
mittivity exhibits the renormalization of the non-Abelian 
coupling coefficients Gauge invariance of the the- 

ory constrains the two renormalizations, as shown in the 
above equations, via the dependence on /eg. Explicitly, 
the two integrals Jq and Ii which define the two permit- 
tivities may be generated from a single function: 



1 d f 

k^duj ^'^i'^)^AH.feqd^P 

ln{V)dpJeqd^p. 



LA 

k^ dk 



(20) 
(21) 



Note that all the "couplings" f^^"^ undergo the same 
renormalization, as required by gauge invariance. 

The resolution of the permittivity into the Abelian and 
the non-Abelian sectors is itself dependent on the gauge 
employed. The above identification is done in the tem- 
poral gauge. However, the results arc themselves gauge 
invariant. We refrain from writing the response functions 
and the counter parts of EqlTH] in an arbitrary gauge as 
it is rather cumbersome, and of no use to us here. 



B. Form of the permittivities 

It is instructive to study the form of the response func- 
tions in the simple case when f^q has the standard Fermi- 
Dirac form. The expression for the Abelian permittivity 
reads 



1 



2TT^Q^T^Nf 
3fc2 



{- 



-In 



cu - 



2} (22) 



and is not different from that of the electrodynamic 
plasma, except for multiplicative color and flavor factors. 
The non-Abelian permittivity is however, novel, and has 
the form 



'^^^^ 9 ZJ^^- 



(23) 



Although the functions are real, their imaginary compo- 
nents can be extracted by employing the standard Lan- 
dau ie prescription. The imaginary components cause 
the Landau damping which we will discuss in detail in a 
separate section. 
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C. Physical significance of ejv 

The significance oi is best seen in the induced non- 
Abelian component of the charge density, the so called 
charge density carried by the field. In the Gauss' equa- 
tion for non-Abclian fields the term pi = f'^^'^A\E'[ 
(where superscript '/' stands for field component) has a 
natural interpretation of the charge density contributed 
by the field. This charge density is not gauge covariant. 
However, the total charge is gauge covariant, provided 
reasonable boundary conditions are imposed, eat may 
now be looked upon as modifying p^ by inducing addi- 
tional charges. To study this we evaluate the modifica- 
tion to pI and obtain 

X I dt'A\{t',r) j dt"E'?{t",r), (24) 

where the second term shows the induced charge density 
of the field. The above equation displays the inherently 
nonlinear, non-Markovian nature of the non-Abelian per- 
mittivity. In the coUisionless limit that we are interested, 
the response is maximally non-Markovian. 

It is pertinent at this stage to mention that the non- 
Abelian response plays an important role involving three 
gluon processes, especially gluonic bremsstrahlung of glu- 
ons, and the analog of Cerenkov radiation of gluons. 
Their phenomenological significance in heavy ion colli- 
sions remains to be investigated. 



D. The gluonic sector 

We now consider the bosonic content of plasma. It dif- 
fers from the matter sector in two respects. Its equilib- 
rium distribution function is given by the Bose-Einstein 
form, and it has a richer structure in the color space. As 
a warm up, we first consider SU(2), which is simpler. 



1. SU{2) gluons 

The color coherent states are no different from the 
more familiar spin coherent states (corresponding to spin 
-1) in the adjoint representation. The density operator 
in the coherent basis has the expansion 



(25) 



where the equilibrium distribution /q corresponds to the 
singlet, /° to the triplet D^, and to the 5-dimensional 
irreducible representation of SU{2). Recall that 
is a completely symmetric and traceless matrix, as is also 
Q°'^ which is given by. 



There are, in all, nine independent functions of position 
and momentum. Of them, the tensor components are 
specific to the gluonic sector and are absent in the quark 
sector. Thus the perturbation around the equilibrium 
configuration has a richer structure for the gluons than 
for the quarks in the color space. We treat /" and 
to be fluctuations around /q. 

We may now repeat the evaluation of the moments in 
the color space as in the quark sector. Since the ten- 
sors are irreducible, the quark sector results hold in the 
gluonic sector as well for the charge density leading es- 
sentially to the same expression for the permittivities 



+ 2}, 



1 - 



1 
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(26) 



(27) 



These differ from the expressions for quarks because of 
the choice of /q. However, unlike in the case of the 
quarks, the above permittivities do not capture com- 
pletely the response of the medium. There are additional 
contributions coming from the tensor fluctuations 
To study them, we evaluate the moment of the transport 
equation wrt Q"^. In the spirit of the earlier calculations, 
we keep only the contribution from /o in the terms in- 
volving the gradient wrt the momentum. In this case, it 
is preferable to evaluate the moment in the phase space, 
i.e., without integration over momentum. Let us denote 
by (!)"■'' {u!, ki, Pi) the Fourier transform of We then 
get 



^ j-abcde^ 



vM"®r"){i^,k,,p,)^0, (28) 



where the symbol ® denote only the convolution of the 
Fourier transforms, without any operation in the color 
space. The fifth rank invariant tensor is given by 

^abcde ^ i(^^o.cdgbe _^ ^bcdgaey ^39) 

The above expression is an eigenvalue equation in the ten- 
sor fluctuations /"'', in the extended phase space. The 
spectrum is, significantly, independent of the equation 
of state, i.e., of /o - within the approximate framework 
employed here. It is determined entirely by the structure 
constants and, of course, the perturbing gauge fields. It 
should be of interest to work out the experimental con- 
sequences, in terms of the propagation of gluons in the 
medium. 

To gain further insight, we integrate Eql^Hl over the 
momentum variables. We obtain 

p'^'icj, k,) - -L-''^^-{A'i ® Jt){uj, h) = 0, (30) 

which establishes a relation betwen the tensor charge 
density 



Q^fj)"^ d^p 
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with the tensor current density 

Observe that Eql3D] bears some resemblence to the conti- 
nuity equation, except that it does not have the standard 
convective term. The absence may be attributed to the 
long wavelength approximation. The corresponding ten- 
sor charge would only be approximately conserved. 



We turn our attention to the tensor fluctuations /"''. 
They satisfy the eigenvalue equations 

^abcd^cd^^^ fc„p,) - 2iL'''''''^%^^A'i ® (j)'^'' = 0, (34) 

where the rank four tensor, specific to SU(3) has the form 



5 

and the rank five tensor has the form 



2. 



The SU{3) gluons 



6 



The generalization to the SU{3) case is not difficult. 
The form of the distribution function is given by 



(31) 



which is formally identical to Eq |25l The transformation 
properties are however different: /g is the singlet, /° the 
octet, and the tensor is the 27-plet of S'C/(3). The 
structure of the tensor basis is more complicated 
than its counterpart in SU{2) [42 |: 



Q2Q°fc = Q^Qb _ hpi-d'-''''QPQ'' - Q2^, (32) 

o 

where d"''^ are the symmetric structure constants of 
SU{3). Note that Q'^^ is symmetric and traceless. In all, 
there are 1-1-8-1-27 parameters that characterize the color 
distribution, as it should be for a 8 x 8 density operator. 
We shall not write the expressions for the Abelian and 
the non- Abelian permittivities since they are identical to 
EqEiand Eqim 

At this juncture we note that there is an interesting 
relation between the permittivities of quarks and gluons 
obeying the ideal equation of state. Consider the sus- 
ceptibilities A{q, g) = — 1 and Af{q, g) = en — 1 for 
the quarks(q) and the gluons(g). It can be easily verified 
that the following relation holds 



Aiq) 



Nf 



A{g), 



mi) = ^AA(g). 



(33) 



It can be shown that it is true for any SU{N) gauge 
group. 



j^abcde ^jacd^be _|_ jbcd^ae _ jmce ^abr ^rdm^ 
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A further integration over momentum yields, in the large 
wavelength limit, 

^abcdpcd^^^ fc,) - 2-L"^^'*^lyl^ = (35) 

which is similar to Eq[30l It has a more complicated 
structure with the occurrence of the symmetric as well 
as the antisymmetric structure constants of the group. 
We study the above equations in a little bit more detail 
in the next section. 

We wish to point out that one needs rather special 
fields to excite the tensor fiuctuations. To demonstrate 
this aspect, let us consider Eg 1551 when the gauge fields 
represent a plane wave, i.e., 

A^UJ, h) - -i^5{uj - u;o)S\h - ho). 

UJ 

For this choice, Eq l35l has only the trivial solution, viz, 
= 0. One may expect that the tensor fluctua- 
tions may be excited provided the field has a quadrupole 
component for SU{2), and its generalization thereof for 
SU{3). The implication of this result to the signatures 
of QGP needs to be explored. 



3. Closure in the color space 

We have seen that the color degree of freedom does not 
possess an infinite hierarchy of moments since its associ- 
ated phase space is compact. For SU{2), the hierarchy 
of the moment equations terminates at the dipole level 
for the matter sector, and at the quadrupole level for the 
gluonic sector, as may be seen below: 



1 cbdfi 



1 



+ iQ^£i^ m—{Ai* 



V 'dp, 



3iV2 1 
Vi v" 



dpi 



= 0, 



(36) 



where Vi , V2 are defined in the Appendix. Similarly, the 
tensor fluctuations satisfy 
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lab 



+ £ 



r 



v ' 



ide 



(37) 



The above equations are written without making the long 
wavelength approximation. This closure as exhibited in 
EqslSni EH is formal since the gauge fields have to be 
determined self consistently with the matter distribution, 



J 



as governed by the Yang-Mills equations. 

The corresponding equations for SU{i) are slightly 
more complicated: 



where the tensor 



Vi V 



dpi 



(38) 



The tensor excitations satisfy the equation 

1 



bed { icd 



rpabcd ^ac^bd _^abe^ecd 

5 



90' 



'V 



dp, 



dpi 



V 



ide 



(39) 



r 



where the two constants Vi and V2 are defined in ap- 
pendix A. We see sagin that the above two equations 
constitute a closed set, in a formal sense. Although their 
full solution in a self consistent manner is an intricate 
and an involved exercise, the clear signal from the above 
equations is that generically, all the fluctuations exist 
and none of the responses can be ignored. It also shows 
that if required, it is possible to improve upon our long 
wavelength approximation in a systematic manner. 



IV. RESULTS AND DISCUSSION 
A. The Abelian response 

Our analysis has made no assumption on the equi- 
librium distribution, save its isotropy in the momentum 
space, and independence of the position coordinates. We 
now consider, purely for purposes of illustration, an ideal 
gas of quarks and gluons. Of course, the Abelian compo- 
nent of the chromoelectric response is no different from 
its electrodynamic counterpart, except for color-spin fac- 
tors. We merely record two of its properties. 

Considering quarks first, we get 



lim e^(w, k) 
lim e^(uj, k) 

k->0 



(40) 
(41) 



written in terms of the standard Debye mass 

ml = AT:^Q^T^Nf/S. 

The corresponding expressions for the gluonic plasma is 
obtained by simply replacing Np by 2 in the expression 
for m|,. It should be borne in mind that since we either 
have a gluonic plasma or a quark gluon plasma, the ex- 
pression for m|) for QGP is obtained by averaging over 
the contributions from the quarks and gluons, weighted 
by the respective relative densities. 

It is commonly accepted that the inverse Debye mass 
yields the screening length for the heavy quark poten- 
tials. This is so for electrodynamic plasmas, governed 
by the Coulomb potential, which goes over to Yukawa. 
The situation is not that simple in the case of quarkonia. 
Since J/5' suppression is one of the strong signatures for 
QGP, it is worth analyzing the import of the Debye mass 
to realistic heavy quark potentials. We take it up in the 
next subsection. 



B. The heavy quark potential 

Consider the Cornell potential given by [i^ 

0(r) = - — + Ar, 
r 

where a and A are phenomenological constants. This 
potential has been studied in depth by Brambilla et al 
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[44| . The potential, has a Coulombic behavior at small 
distances and a linearly confining form at large distances. 
Note that the parameter A has mass dimension 2 (a has 
mass dimension 0). The medium alters the above poten- 
tial as given by the modified expression for the potential 
(in the Fourier space) 4's{k) = (j>{k)/e{k). 4)s{k) 

is obtained to be 



given by 



<Ps{k) 



A 



7rfc2 



-D 



277 +m|,) 



(42) 



The inverse Fourier transform of the above expression 
yields the modified potential, as a function of distance, 
as shown below: 



Ur) = (^ 



a) 



exp{—mDr) 2A 



(43) 



The above expression merits a closer scrutiny. Interest- 
ingly, the Coulomb term goes over to the short range 
Yukawa form while the linear confining term goes over to 
a sum of Yukawa and the long range Coulomb potential. 
The Coulombic tail, which is relevant at large distances, 
supports - as is universally known - an infinite num- 
ber of bound states. The effective Coulombic charge is 
given by ^4-; this suggests that the role played by the 

Debye mass cannot be naively reduced to the notion of 
a screening length. Consequently, the relevant parame- 
ter to study the deconfinement is, not the inverse Debye 
mass, but the dissociation energy as a function of the pa- 
rameter -^4-. Note that it is also the relevant parameter 

at high temperature. 

Thus in the large distance/high temperature limit, it is 
not difficult to estimate the dissociation energy £d- We 
obtain 

£D-m,A^/mj,^-^-I^. (44) 



(47r 



It is further possible to estimate the temperature To 
at which the dissociation takes place. Employing the 
equipartition theorem, we obtain 
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(47r3)2 NfQ^ ■ ^"^^^ 

We note that for a purely gluonic plasma, the factor Np 
gets replaced by 2. 

In short, a rather straightforward and a simple analysis 
has allowed us to extract the precise physical significance 
of the Debye mass for the dissociation, or equivalently, 
the suppression of the charmonium states in QGP. It also 
demonstrates that care should be exercised in interpret- 
ing the inverse mass as a screening length. The physics 
is dictated by a combination of two scales, A and m|,, as 
obtained in EaH51 



C. The non-Abelian response 

We turn our attention to the non-Abelian permittivity. 
In the long wavelength limit, it has a simple expression 



Cat 



1 



It is interesting to note that the above expression is al- 
most identical to the expression for the Abelian permit- 
tivity at fc = 0. In fact, the quantity ( in eA gets 
simply replaced by ucu' in e^r. The non-Abelian response 
vanishes at large a;,w', and the response gets stronger 
with temperature. Of course, it vanishes when T = 0, 
but it cannot be taken seriously since the interaction ef- 
fects would dominate at small temperatures. As men- 
tioned, we may expect the non-Abelian response to play 
significant role in processes involving three gluon vertices. 



V. LANDAU DAMPING IN QGP 

Landau damping in QGP has been studied earlier by 
Heinz and Siemens [45[ and Murtaza, Khattak and Shah 
[46| . We obtain here an explicit expression for Landau 
damping in QGP, arising especially from the non-Abelian 
response. A naive extension of the standard formula will 
not do since we have additional permittivities. Further 
the Abelian and non-Abelian components of the permit- 
tivities are not independent because of gauge invariance. 
We address the problem ab initio. 

The interaction Hamiltonian is given by 



(46) 

The energy dissipation per unit time per unit volume has 
the standard expression 

Recall that we are in the temporal gauge. Employing 
Ea ll21 which we reproduce below, 

1 dfo 



it is straight forward to obtain the contribution to the 
energy density due to the non-Abelian response. Finally, 
adding the standard Abelian contribution, the total en- 
ergy density per unit time of an electric field is obtained 
to be 

U = uj{l + Q^Ioiu;,k)}E,-E, 

+Q^h{uj,k)\k=oE,-E,. (47) 

In Eg 1471 the contribution from imaginary term gives the 
rate of the energy dissipation rate per unit volume, Sdis- 
The contributions from Iq and Ji are respectively from 
the Abelian and the non-Abelian permittivities. It is 
noteworthy that in the coUisionless limit that we are in- 
terested in, only the Abelian component contributes to 
the damping, for, Eq[47| leads to the explicit form 



lld^s = T^mll^eik + ^}E, ■ E,. 

The non-Abelian contribution is saturated at lu 
is hence not an observable. 



(48) 
and 
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VI. CONCLUSION 

In conclusion, we have developed a systematic for- 
mulation for determining the chromoelectric response of 
QGP within the classical framework. We have shown the 
emergence of a non-Abelian component of the response, 
characterized by new permittivity cn. ejq is, unlike 
its Abelian component e^, nonlocal and non-Mar kovian. 
Furthermore we have also derived explicitly the response 
functions of gluons which are shown to be richer than 
that of the quarks. The precise role of the Debye mass in 
the so called screening of the heavy quark potential is ex- 
tracted. Nevertheless, we point out that the applications 
discussed in this paper are only indicative in nature since 
we have ignored the strongly interacting nature of QGP; 
our expressions for the equilibrium distribution functions 
are highly approximate. An application of the analysis 
made in this paper, with realistic equations of state as 
applied to QGP will be taken up separately. 



VII. APPENDIX 
A. Color space integrals 

The Haar measure for the SU(3) group is given by 



[13, where ~ N and = in the adjoint representa- 
tion of SU(N) group [30|. Some useful color integrations 
in SU(3) group space are given by 



Q^Q'''' dQ = 0; I Q^Q' dQ = ^Q^i; 



where Vi is the volume of color space and V2 is a constant. 
Q"** has been defined as (see Ea l32)) 



r'' = Q^Q^ - -dP«'^d™''QPQ'? - —Q^. 
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